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The Pomeron-gluon-gluon interaction is considered in the QCD- 
CN . based model for the charmed baryon production in the process 

Pomeron + p — > + X. The polarization of the produced 
heavy quark is induced effectively through the non-perturbative 
long-range interaction with the gluon field of the type [a ■ rot A] . 
The ^-dependence of A+ polarization, P A +(xp,pp), has been 
studied. Its absolute value depends on the model parameter a 
qq ! and it appears to be sizeable in the wide range of a values: when 

a ranges from 0.1 to 1.0, the polarization P a +(xf,Pt) varies from 
^ -0.2 to -0.5 at xp ~ 0.5 and pp in the interval l-i-2 GeV/c. 

o ■ Introduction. Polarization properties of heavy quarks have been studied 

for a long period, both experimentally and theoretically. A major piece of 
knowledge came from the experimental study of strange A baryons. The 
polarization of the final-state A's at large xp is negative ant it reaches 20- 
30%. The most recent and full review on the polarization of the strange 
■ baryons can be found in K. Heller's report at SPIN96 JTJ. 

So far, the experimental data on the charmed A+ polarization have been 
restricted to the output of only one experiment performed in IHEP, Protvino 
at the neutron beam in the reaction n + C12 — > A+ + X 0] . The measured 
upper limit of the absolute value of A+ polarization was equal to 0.5 ±0.2 . 
The momentum of the neutron beam ranged from 40 to 70 GeV/c. 

Perhaps soon new data on the A+ polarization will come from the E-831 
experiment at Fermilab (FOCUS) in which the polarization measurements 
can be performed in the reaction 7 + p — > A+ + X at 300 GeV/c. 

The forthcoming collider experiments at RHIC and LHC will make it pos- 
sible to study the A+ production and its spin properties in various kinematic 
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regions, in particular, in the diffraction region. Theoretically the spin phe- 
nomena in diffractive processes are under intensive study. Three refs. in || 
represent different aspects of the spin dependence of diffractive interactions. 

Unlike the strange A, the theory predictions for the A+ polarization P are 
not numerous. The large values up to 100% were predicted for the absolute 
value of P A + in an early paper The recent estimates were made in ref. 
H for the intrinsic charmed sea contribution to the A+ production. Under 
various model assumptions, P A + ranges from -0.02 to -0.11 at large xf in 
the recombination approach. 

Note also that the spin of A+(2285) has not firmly been established so 
far. In the following we shall assume that the spin is equal to 1/2 as it is 
accepted by the PDG §. 

In this paper we are interested in the A+ polarization in the diffractive 
production process 

p + p^A+ + X (1) 

at high energies (fig. 1). In the upper block the cc pair is produced and c 
quark is polarized in the final state. The model for the heavy quark polar- 
ization depends on the kinematics in the upper block and does not promptly 
depend on the diffractive production of the upper block. The advantage of 
the diffractive cc production in (|l|) displayed in fig. 1 is the simplicity of the 
lowest-order model with the dominance of the Pomeron-gluon-gluon coupling. 
The more general consideration would require an additional determination of, 
at least, the ratio of the coupling constants Pomeron-gluon-gluon/Pomeron- 
quark-quark. 

Formulation of the model. We study the inclusive production of cc 
pair by the initial Pomeron in the process 

p + Pomeron -> A+ + X (2) 

(fig. 2) which is relevant to the inclusive reaction (|IJ). In this reaction the 
incoming proton energy \fsj2 and the Pomeron 4-momentum q are those 
of reaction (|T|). q is fixed at an arbitrary small value. The results of the 
calculations of the polarization P a +(xf,pt) do not depend essentially on q 
in the diffractive region. 

The underlying parton subprocess is shown in fig. 2 as the lowest-order 
heavy quark pair production by the gluon-Pomeron coupling: 

g + Pomeron — * c + c. (3) 
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We neglect the charmed sea, so the leading order process c + Pomeron — > 
c + g does not contribute to the A+ polarization. Also we neglect parton 
transverse momenta. The A+ transverse momenta in are taken in the 
region pp = 1 -j- 2 GeV/c. 

Another assumption is the dominance of the Pomeron-gluon-gluon cou- 
pling constant gp gg compared to the Pomeron-quark-quark coupling gp qq , so 
that gp q q/gpg g <C 1. Hence, at this level we exclude from the consideration 
the next-order subprocess q + Pomeron — > q + cc in which a valence quark 
q is excited by the Pomeron and it emits the cc pair. 

The produced c quark recombines inclusively with the spectator ud di- 
quark system thus forming the final A+. 

The A+ polarization originates from the polarization of the produced c 
quark and is expressed in terms of the quark spin density matrix pg A by the 
following formula for the charmed baryon spin density matrix in process (Q): 

p AA ' = J dx ud dx g dz f g {x g , Q 2 ) fud(x ud , Q 2 ) £q|| p Q x '(p Q , A) 

xR(x F ,p T ,x Q ,x ud ;z). (4) 

Here f g and f ud are gluon and diquark distribution functions with x g and x u d 
standing for the momentum fractions of the gluon and diquark, respectively. 
Q 2 is the QCD scale chosen in the form Q 2 = m 2 ^ + p\. The quantities 
a, Pq, Eq refer to the underlying subprocess (|). The hadronization process 
in ([|) is described by the inclusive recombination function R depending on 
the c quark, diquark and A+ momenta with z being the energy fraction of 
(ud,c) system carried away by A+. Using (^) the A+ polarization is written 

as P A+ = Tr ( (J P)/ Tr (p)- 

Initially a model for the final-state quark polarization expressed in terms 
of pQ X (pq,A) was formulated in 0. In this paper we explore a corrected 
version of that model applied to the charmed baryon production. In the 
model the produced heavy quark scatters off the external gluon field A and 
the quark spin sq may couple to the field, so that the correlation sq ■ rotA 
occurs. 

To estimate this correlation, a non-trivial model for the external gluon 
field A itself is needed. In case of the simplest color field, the spin flip 
properties and polarization of a scatttered massive quark were considered, 
for example, in the early paper || by J. Szwed, where the external gluon field 
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$ a (g) = Airgl a /(f was used in the perturbative second-order calculations. 
Note that the massive quark polarization obtained in || was negative. 

An alternative way is to parametrize the heavy quark spin density matrix 
Pq~ globally rather than to introduce an explicit model for A and for the 
quark interaction with this field. 

The parametrization of pq is based on the following considerations. By 

definition Pq X = L Xfl L x M where L Xfl is the quark transition amplitude 
in the gluon field A. The following traces relate to the basic formula (H): 

Tr(o-p Q ) = -2ImpJ- = 41m (L ++ L + ~) and Tt(pq) = 2(|L++| 2 + \L+-\ 2 ). 
Under the space parity conservation the relations L ++ = L and L + ~ = 
-L-+ hold. 

For the parametrization it is convenient to include Tt(pq) into the recom- 
bination (hadronization) function R in (Q). Then the trace of the convolution 
of (^) with the Pauli matrix a will contain in the integrand the expression 
Pq = Tr(crpQ)/Tr(pg) which is exactly the heavy quark polarization (see 
below). 

After these generalities the time has come to parametrize the ratio Pq. 
The normalized spin-flip transition amplitude squared |L H | 2 has the sense 
of probability to make quark spin-flip in the external field. It seems natural 
to think that it increases with the quark mean range. Henceforth, it grows 
with increasing quark energy, Eq, since at hadronization the mean range of a 
fast quark is proportional to Eq/M 2 where M _1 is of order of the confinement 
radius. Thus, we accept the equality |L + ~| 2 /(|L ++ | 2 + \L + ~\ 2 ) = a£Q / EQ ax 
with the supression parameter 0< a <1. The normalization EQ ax depends, 
of course, on the kinematics, however, for the simplicity we take it equal to 

Taking these straightforward considerations into account the parametriza- 
tion has been made as follows: 

2Im(L++ L 



P Q = | L ++|2 + | L +-|2 =2^/aa; e (l-aa; e )sinA0, (5) 

where x £ = Eq/eq"^, and A0 is the phase difference of the amplitude L ++ 
and L + ~. 

Phase difference. In principle the phase difference A0 depends on 
arguments of the amplitudes L xx , i.e. on Eq. The simple linear dependence 
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A0 = tiEq + r 2 was used in || to describe the experimental data on the one- 
spin asymmetry in the reactions p-p + p — > tt^ + X at 13 and 18 GeV/c and 
in the reaction tt~ + p^ — > 7r° + X at 40 GeV/c. Three sets of experimental 
data were enough to determine the parameters r\ and r 2 . 

In this paper we prefer to imply A<p value based on the Regge-type con- 
siderations. The Regge properties of the reaction (|2]) are determined by D* 
and D** exchange trajectories (see insertion in fig. 1). And in the Regge 
approach the A+ polarization in reaction (0) is due to the interference of 
the exchange diagrams with these trajectories J3J. Thus the difference of 
the Regge phases, A0 = <pp>* — <fi D ** occurs naturally. It depends on t, 
the squared momentum transfer from the initial proton to the final charmed 
baryon. Since we consider A+ production at large xp and relatively small p?, 
we neglect the small t contribution in A</> (in || the t-dependence dropped 
down due to the hypothesis of the weak exchange degeneration). So, we ob- 
tain Acj) = ±(j* — j**)n/2 = ±7r/2. Here the starred j's correspond to the 
spins of the exchange trajectories. The sign ± reflects the non-uniqueness 
of this procedure and it can be chosen by comparing the results with the 
strange A polarization at large xp. 

Parton distributions. We used two sets of the leading-order parton 
distribution functions (pdf's): CTEQ4L and GRVLO fpLOf] - 111 the kinematic 
region of interest, the gluon pdf's f g (x, Q 2 ) from these sets are close at x < 
0.3 -r 0.4 and they differ at larger x, whereas the pdf's of the light quarks u 
and d are similar. 

For our purpose, the 'diquark' pdf has been constructed of u and d quarks 
very simply: 

J y x-y 

The arguments in the denominators means that the pdf's output from the 



packages in |10| are of the conventional type x times pdf. The diquark distri- 
bution is shown in fig. 3. 

Inclusive recombination. The hadronization of the produced charmed 
quark is assumed via the inclusive recombination with the spectator valence 
ud system. This subprocess is described by the inclusive recombination func- 
tion 

R(x F ,p T ,x Q ,x ud ;z) = const ^— 8(z - 1) 

•Esq ~r -^ud % 
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where xe is the reduced A+ energy: xe = 2E/y/s, and z = xe/(x Eq + x u d). 
This expression incorporating the transverse momentum is a generalization 
of the conventionally used recombination functions. The normalization pa- 
rameter const has not been specified since it drops down in the formula for 
polarization P A +. 

Partonic subprocess. The differential cross section 6Q(l 3 a /dpg has 
been calculated using the matrix element 

M = g Pgg g s U P T^e Ql ^ Q . 

Here gp gg and g s are the Pomeron-gluon-gluon and the strong coupling con- 
stants, respectively, is the gluon polarization 4-vector, uq and vq are the 
quark and antiquark spinors, respectively, T£ is the color factor and s is 
the incoming energy squared. (pp is the Pomeron wave function which is 
normalized as 4>p(pp=l. 

Kinematics. We are interested in xp dependence of the A+ polarization 
P a +(x p ,pt) in the process (0) with given p T and at a fixed value of q. The 
kinematic region 0.2< xp <0.9 and 1< pp <2 GeV/c was studied. No visible 
Pp dependence of P a +(xf,Pt) was obtained. 

The basic relations for the longitudinal and transverse momenta are the 
following: xp = x u d + Xq, pp = ppq. This should be added by the 4- 
momentum conservation in the partonic subprocess. Thus, the integration 
limits in x u d and X g CcLll be formally determined. However the lower limit in 
x u d, the momentum fraction of the spectator valence diquark would be too 
small from the physical point of view. So the lower limit of x u d was chosen to 
be not less than 0.15. Variations around this value do not cause the strong 
change of the resulting polarization. 

The lower limit in integrating over x g , the gluon momentum fraction, is 

defined by x g nm = (2/ ^/s){\Jpp + rriq + sq — go)- The value of the upper 
limit has been chosen as x 1 g nax = 0.95. 

An important note. Here is an appropriate place to make a note on the 
seemingly numerous assumptions before starting to calculate the polarization 
integrals. The final-state baryon polarization at large xf is due entirely to the 
QCD twist-3 contributions. For the twist-3 studies in hadron interactions, 
there is no commonly used model. At first glance, there is too much freedom 
in the above considerations of the quark polarization Pq and/or the inclusive 
recombination function R. It seems to be a disadvantage of the approach. 
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However the hadron (A+ in this paper) polarization is a ratio of two integrals 
which include the same functions in the nominator and denominator, except 
Pq presented in the nominator only. And all the functions are rather smooth. 
From this it follows that the ratio of the integrals depends weakly on the 
reasonable variations of the separate functions in the integrands. It means, 
for example, that our results must not change strongly if we use the second 
order diagrams in eQda/dp 3 instead of the leading order calculations. The 
same is true for different forms of the inclusive recombination function R. 
The shape and size of the A+ polarization are mostly determined by the 
parametrization of Pq in (|5|). 

Results of calculations. As is seen from (|j), the sign of the A+ polar- 
ization is determined by the sign of Tr(crpQ), that is the sign of the heavy 
quark polarization Pq. The presented model for Pq (^) has the sign ambigu- 
ity mentioned above in the paragraph Phase difference. Since the model 
makes no difference between the flavors of the massive quarks, we resolve the 
sign ambiguity in favor of the strange s quark in order to obtain the negative 
A polarization at large Xp in similar calculations. So the phase difference 
A(p in d5|) has been chosen equal to -7r/ 2. It is worth to note that this choice 
supported by the experimental data for the strange A polarization coincides 
with the sign of Pq in || obtained in perturbative calculations. 

The only undetermined quantity is still the suppression factor a in formula 
(|5]). The numerical calculations were made for two far-distant values a=0.1 
(small polarization) and a=l (largest polarization) using two sets of the 



leading-order parton distributions CTEQ4L and GRVLO P0[ . The results 
for xp dependence of the A+ polarization P a +(xf,Pt) at px= 1.5 GeV/c are 
shown in fig. 4. At xf=0.5 the small and the largest polarizations differ in 
magnitude by a factor of ~2.5. With growing xp from 0.2 to 0.9, both the 
small and the largest polarizations increase by a factor of 2. 

Acknowledgement. We would like to thank Hung-Liang Lai, one of the 
authors of CTEQ4 pdf package fr0"|, for the instructive help in using CTEQ. 
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Fig. 1. (a) - Diffractive production of with cc pair in the upper 
block in reaction (1); (b) - Regge exchanges in the upper block. 




Fig. 2. The leading-order parton subrocess in reaction p + 
Pomeron — ► A+ + X. 
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Fig. 3. Diquark distribution function xf u d(x, Q 2 ) at Q 2 
with rriQ — 1.5 GeV/c and pr = 2GeV/c. The solid (dashed) line 
corresponds to CTEQ4L(GRVL0) parton distribution functions. 
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Fig. 4. xf dependence of the polarization in the process p + 
Pomeron — > + X at pr = 2GeV/c. The solid (dashed) line 
corresponds to CTEQ4L(GRVL0) parton distribution functions. 
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